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Method and associated device for generating random 
numbers in a given range 

The invention concerns a method of obtaining a 
random number between A and B from a generator 
5 producing random numbers lying between 0 and W-l, with 
N the size of the numbers produced by the generator, W- 
1 the maximum value taken by the random numbers 
produced, with for example W = 2 N and A, B any integer 
numbers, less than or greater than the number W. 

10 Such a situation occurs for example in an 

electronic component adapted to perform cryptographic 
calculations and comprising an N-bit random number 
generator, for example N = 8. The random numbers that 
it can produce are thus between 0 and W-l = 255, whilst 

15 it would be desirable to have random numbers between 
for example 0 and 100 or between 300 and 10000. It 
should be noted that it suffices to determine numbers 
between 0 and 9 70 0 and then to add 3 00 to the number 
obtained in order finally to obtain a number between 

20 300 and 10000. 



Such a situation is found in practice in the 
majority of cryptographic applications, for example the 
DSA signature, the El Gamal signature or enciphering, 
the development of countermeasures against various 
attacks, etc. 

Several methods are already known for producing 
random numbers R between 0 and K from numbers between 0 
and W-l. These methods are in general implemented by 
software means used to control on the one hand a 
hardware generator that produces random numbers of size 
N and on the other hand calculation means performing in 
particular multiplication, addition, etc operations. 

A first known method comprises the following 

steps : 

a) determining the smallest integer number p such 
that K < WP - 1, 

b) producing p random numbers S 0/ S, S p _! and 

forming the variable S = ^S i *W i 

i=0 

c) if S > K, then returning to step b) , otherwise 
putting R = S 

R is the random number sought, between 0 and K. 

The equation S = ^S i *W i is a representation of the 

i=0 

variable S decomposed/recomposed in base (W p_1 , . w 1 , 
w °) • It would also be possible to note S 

Sp-iSp_2 . - . SiSo / a notation commonly used. 

A second known method comprises the following 

steps : 



a) determining the smallest integer number p such 
that K < WP - 1, 

b) producing p random numbers S 0/ S, . . . , S p _i and 

forming the variable T = 2^S i *W and S = T + Sp-i*!^ 15 " 1 

i=0 



c) if S > K, putting R = T otherwise putting R = 



A third known method comprises the following 

steps : 

a) determining the smallest integer p such that K 
< WP - 1, 

b) producing p random numbers S 0 , S, . .., S p -i and 
forming the variable S = ^ j S i *W l 



c) putting R = S mod(K+l), that is to say the 
remainder of the whole-number division of S by K+l, 
also referred to as modular reduction of S by K+l. 

These three methods can be summarised by the 
following steps: 

a) producing p random numbers S 0 , S, . S p -i, 
being the smallest integer number such that K <W P - 1 

and forming the variable S = ^TS i *W i 



b) determining the random number R from the 
variable S . 

According to circumstances, during step b, R is 
obtained from S by repeating step b (first method) , 
taking account or not of the additional random number 
Sp-i (second method) or performing a modular reduction 
(third method) . 
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It should be noted that, in the three methods, if 
a number between A and K+A is required, it suffices to 
add A to the number R obtained lying between 0 and K. 

The main drawback of the first method is a 
5 particularly long and especially unpredictable 
calculation time: the step of producing the p random 
numbers may be repeated numerous times without it being 
possible to predict at the start the number of 
repetitions of this step. 

10 The second and third methods have the main 

drawback of producing random numbers exhibiting a bias: 
amongst the numbers R produced in the range [0, K] , 
certain values are more probable than others. In other 
words, the numbers R produced are not perfectly random 

15 (non-uniform distribution) . This bias may have 

significant consequences on the security of the 
cryptographic systems liable to implement these methods. 
The security of cryptographic systems assumes in fact 
that the random numbers that they use are uniformly 

2 0 distributed (or at least close to a uniform distribution) 
in the range [0, K] or [A, K+A] wished for. 

Finally, the three methods are slow overall 
because they implement operations on large numbers, of 
size N (in the sense of the number of bits) greater than 

25 the size of the circuits used for the implementation. 

This is because the number K in particular is any number 
and can be greater than W and therefore of size greater 
than N. The variable S can also be of large size. 
However, the implementation of operations on large 
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numbers requires the implementation of complex methods 
expensive in terms of calculation time. 

An essential object of the invention is to 
propose a method of constructing a random number R that 
5 is particularly rapid. 

Thus the invention proposes a cryptographic 
method during which use is made of a random number 
generator producing random numbers Si of size N fixed 
between 0 and W-l, with for example but not necessarily 
10 W = 2 N , in order to produce a random number R between 0 
and a predefined limiter K. 

The essential steps of a method according to the 
invention are as follows: 

E31: a random variable Si between 0 and W-l is 
15 produced, 

E32 : if the random variable Si is strictly less 
than a coefficient Ki of the limiter K in base W, then 
the coefficient Ri of rank i of the random number R is 
equal to the random variable S ± and then, for any rank J 
2 0 less than i, a random variable Sj between 0 and W-l is 
produced and Rj = Sj , 

E33: otherwise, if the said random variable is 
greater than the coefficient Ki of rank i of the limiter 
K in base W, then the said coefficient R± is determined 

2 5 from the random variable Si of rank i according to a 

predetermined function, and then the coefficient Ri- X is 
determined for the random number R of rank i-1 that is 
immediately lower by repeating steps E31 to E33 . 

Thus, in a method according to the invention, the 

3 0 coefficients Ri of the random number R required are 
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sought one by one, commencing with the most significant 
coefficient R p -i . The physical generator of random 
numbers used thus produces random variables Si one by 
one, one variable at each iteration. 
5 In addition, the method is rapid since step E33 is 

executed a small number of times. This is because, as 
soon as one of the variables Si produced by the physical 
generator is less than the associated coefficient Ki of 
the limiter K, the method no longer requires the 

10 processing of the variables Sj of rank less than i: thus 
a small number of coefficients of the number R, the most 
significant, are calculated the most often. 

Finally, compared with the known methods, a 
method according to the invention has the advantage of 

15 working on numbers of no more than N bits, N being the 
size of the registers and other calculation circuits of 
the devices used for implementation. For example, if W 
is equal to S N , the coefficients Ki resulting from the 
decomposition of K in base (W p_1 , ... W 1 , W°) are 

20 necessarily less than W and therefore with a size of no 
more than N bits. Likewise, the random variables Si 
produced by the physical random number generator are 
also of N bits. 

By adding to the essential steps an 

25 initialisation step and a step of recombination of the 
random number R, there are obtained: 

El: the limiter K is decomposed in base (W p_1 , W p " 2 

W°) (K = f^K i *W i or K = K p " 2 ... K 1 ^) , i being a 

i=0 

loop index, K± being a coefficient of the limiter K of 
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rank i between 0 and W-l and p being the degree of the 
limiter K, 

E2 : a Boolean variable f is initialised to TRUE, 
E3 : the following operations are performed, in a 
5 loop indexed by i, i being an integer varying between 
p-1 and 0: 

E31: a random variable Si between 0 and W0-1 

is produced, 

E32 : if the random variable Si is strictly 

10 less than the coefficient K± of rank i, then the 

Boolean variable f is set to FALSE, 

E33_l: if the random variable Si is strictly 

greater than the coefficient Ki of rank i and the 

Boolean variable f is TRUE, then the coefficient Ri 

15 of rank i is determined from the random variable Si 

of rank i according to a predefined function, 

E33_2: otherwise Ri = Si 

E34 : the loop indexed i is decremented, 

E4 : the random number R is determined by 

20 recombination of the random coefficients Ri in base 
P -\ 

W{R = Ya R * wi or R^.^RV). 

i=0 

In concrete terms, as soon as the Boolean 
variable f is positioned at FALSE, it remains at this 
value since provision is not made for repositioning it 
2 5 at the value TRUE, except when E2 of the method is 
initialised. Step E32 is executed only if the variable 
f is TRUE; thus, as soon as the variable f is 
positioned at the value FALSE, step E33_l is no longer 
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executed and the method according to the invention ends 
rapidly. 

A second objective of the invention is to propose 
a method of constructing random numbers whose 
distribution is uniform or can be made as close as 
desired to a uniform distribution. This objective is 
achieved by choosing a suitable function for the 
determination of the coefficient R ± from the random 
variable Si. 

According to a first embodiment of the method 
according to the invention, in order to determine the 
coefficient R ± of rank i from the random variable Si of 
rank i (step E33_l) , the following substeps are 
performed: 

E33__ll: if the random variable Si is strictly 
greater than the coefficient K t of the limiter K, then a 
new random variable Si is produced, 

E33_12: step E33_ll is repeated until the random 
variable Si is less than the coefficient K± of the 
limiter K, and then the coefficient R ± is equalised to 
the random variable Si. 

In such an embodiment, all the coefficients R± 
obtained are numbers directly produced by the hardware 
random number generator; and these coefficients are 
therefore perfect and the number R which results 
therefrom is also perfect. In other words the 

distribution obtained of the numbers R is uniform in 
the range [0, K] . 

According to a second embodiment, during step E33 
the coefficient R± of rank i is chosen so as to be equal 



to part of the random variable Si, a part less than the 
coefficient Ki. The said part corresponding in one 
example to a limited number of bits of the variable Si. 

According to a third embodiment, during step E33 
the random variable Si is reduced modulo Ki + 1, the 
results of the reduction being the coefficient Ri 
sought . 

These latter two embodiments are rapid compared 
with the known methods, essentially because the work is 
done on small numbers. The distributions of random 
numbers obtained are however not uniform: the simple 
fact of truncating the variable Si or performing a 
reduction modulo Ki+1 necessarily introduces a bias. 
However, this bias is less compared with the methods of 
the prior art . 

Moreover, it is possible to reduce the bias of 
the methods according to the second and third 
embodiments proposed, as will be seen below. 

In a method according to the invention as 
described above, a random number R is constructed less 
than K from variables S± of size N produced by a 
perfectly random physical generator. The number R 
obtained is biased, but the bias is small compared with 
a known method. 

For this, in the second or third embodiment, a 
coefficient R ± < K ± is constructed, in particular during 
step E33__l, from variables Si of size N. In order to 
reduce the bias introduced on the coefficient R ±/ it is 
proposed to construct it using the same steps El to E3 
as for constructing the number R. In a sense, two 



similar methods are "interleaved" . This makes it 
possible to reduce further the size of the numbers on 
which the work is carried out, and consequently to 
reduce further the bias on the coefficient of R, and on 
the final number R. 

In concrete terms, in order to determine the 
coefficient Ri of rank i from the random variable Si of 
rank i (step E33_l) , steps El to E4 are executed using 
a base (P q_1 , . .., p°) as the calculation base, p being 
an integer number strictly less than W and q being the 
degree of Ki in base p . 

Step E33 is thus broken down into the following 
substeps : 

E33_41: the coefficient Ki of rank i of the 

limiter K in base (P q -\ . .., p°) ( K x = ]T(/^) y * J3 J or Ki = 

(K±) q -i. . . (Ki) i (Ki) o) , j being a loop index, (Ki)j being a 
number between 0 and p-1 and q being a degree of the 
coefficient Ki, is decomposed, 

E33_42: a second Boolean variable g is 
initialised to TRUE, 

E33_43 : the following operations are performed, 
in a loop indexed by j varying between q-1 and 0: 

E33_431: a random variable (Si)j between 0 
and p-1 is produced, 

E33_432: if the random variable (S ± )j is 
strictly less than the coefficient (Ki)j, then the 
second Boolean variable g is set to FALSE, 
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E33_4331: if the random variable (Si)j is 
strictly greater than the coefficient (Ki)j and the 
second Boolean variable g is TRUE, then a 
coefficient (R±)j is determined from the random 
5 variable (Si)j according to a predefined function, 

E33_4332: otherwise, (R±)j = (Si)j 
E33_434: the loop indexed j is decremented, 
E33_44: the random number Ri is determined by 
recombination of the random coefficients (R±) j in base P 

10 or R ± = (Ri) q -i. . . (Ri)i(Ri)o) - 

As has just been seen above, by "interleaving" 
two methods, the bias of the random numbers R produced 
by the global method is reduced, whilst preserving a 
rapid global method. It is of course possible to 
15 imagine "interleaving" more than two methods, for 
example three or four, by decomposing, in step E33_43, 
the numbers in base y < p, and decomposing step E33_43 
in a succession of steps similar to steps E33__41 to 
E3 3_4 3 . 

2 0 In general terms, the more methods are 

"interleaved" , the smaller the numbers on which the 
work is carried out: the duration of each step 
decreases and the bias of the numbers produced by the 
global method also decreases. 

25 Another object of the invention is an electronic 

component adapted for implementing the method as 
described above. Such a component comprises in 

particular a generator producing random numbers of size 



N, and calculation circuits for performing operations 
on numbers of no more than N bits. 

According to the embodiment of the method to be 
implemented, the calculation circuits are adapted to 
perform operations of comparing two numbers, number 
truncation and modular reduction. 

The random number generator and the calculation 
circuits are preferably controlled by a software means 
stored in a memory of the component provided for this 
purpose . 

The invention also concerns a chip card 
comprising an electronic component as described above. 



